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ABSTRACT 

We have observed molecular-line transitions between 84 and 147 GHz in 24 dense 
southern molecular clouds with declinations south of -30 degrees using the Mopra and 
SEST radio telescopes. The results of observations of 2-mm and 3-mm transitions of 
the molecules CO, CS, C2H, CH3OH, HCN, HNC, HCO+, HC3N, OCS, HNCO and 
SO and several of their isotopomers are discussed in this paper, including ratios of 
the emission from the ^^c, ^^C and ^^O isotopomers of CO, CS, HCO+, HCN and 
HNC. This paper lists the calibrated observational data for each molecular cloud 
and discusses some general properties of this sample. The calibration process for the 
Mopra antenna of the ATCA is described in an appendix. 
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1 INTRODUCTION 



Molecular-line emission provides a powerful method for in- 
vestigating the events associated with star formation. Molec- 
ular emission can be used to find the temperature, density 
and kinematics of gas from which stars form, while emission 
from different molecules traces gas of differing physical con- 
ditions. The theoretical and observational underpinning for 
this type of study has been the work of many investigators 
over a thirty year period; a good review can be found in van 
Dishoeck & Blake (1998). 

The southern part of the Galactic plane includes the in- 
nermost parts of the galaxy as well as major dense molecular 
clouds and star-forming regions, many of which are difficult 
to observe with northern hemisphere telescopes due to their 
southern declinations. Consequently molecular line emis- 
sion, molecular abundances and physical conditions have not 
been extensively investigated in clouds at declinations south 
of -30 degrees. 

In the radio spectrum, molecular lines occur primarily 
at millimetre and submillimetre wavelengths. The commis- 
sioning of a millimetre wave receiver on the Mopra ^ antenna 
of the Australia Telescope in 1994 provided the opportunity 
to observe molecular transitions at 3 mm with high veloc- 



ity resolution, making it possible to produce a statistical 
database of molecular spectral-lines and molecular abun- 
dances in southern molecular clouds. 

In the study reported here, a sample of 24 molecu- 
lar clouds with declinations south of -30 degrees has been 
observed with the Mopra Telescope in the 3-mm transi- 
tions of the molecules CO, CS, C2H, CH3OH, HCN, HNC, 
HCO+, HC3N, OCS, HNCO and SO and several of their iso- 
topomers. In addition, complementary observations of 2-mm 
and 3-mm transitions have been obtained with the Swedish- 
ESO Submillimetre Telescope (SEST) ^ 



2 OBSERVATIONS AND EQUIPMENT 
2.1 Equipment 

2.1.1 The Mopra Telescope 

The 22-m 'Mopra' antenna of the Australia Telescope (op- 
erated by the Australia Telescope National Facility, CSIRO 
and the University of New south Wales) is situated near 
Coonabarabran in New South Wales, Australia. The tele- 
scope is at latitude 31 degrees south and has an elevation 



^ The Mopra telescope is part of the Australia Telescope which 

is funded by the Commonwealth of Australia for operation as a ^ The SEST is operated jointly by the Swedish Natural Science 

National Facility managed by CSIRO. Research Council and the European Southern Observatory. 
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of 850 metres above sea level. In 1994 the telescope was fit- 
ted with a millimetre-wave receiver allowing it to operate 
at frequencies between 85 and 116 GHz. The antenna has a 
'shaped' parabolic surface area to increase its forward gain 
and until 1999 it had a solid reflecting surface extending 
over a diameter of only 15 metres. Thus at millimetre wave- 
lengths the reflecting surface had an effective diameter of 15 
m, resulting in a half-power-beam-width that varied from 
47 arcsec to 35 arcsec across the range 85 - 116 GHz. The 
antenna was resurfaced in 1999, extending the effective 
millimetre-wave diameter to 22 metres, giving a half-power 
beamwidth varying from 35 to 24 arcsec over the same range 
of frequencies. 

The Mopra mm-wave receiving system includes a 
cryogenically-cooled SIS (Superconductor Insulator Super- 
conductor) amplifier that can be tuned to either single or 
double side-band mode. For single sideband operation, the 
response of the unwanted sideband can be decreased to at 
least 10 dB below the other sideband response using a cal- 
ibration signal injected into both sidebands. A polarization 
splitter is used to separate the incoming signal into two chan- 
nels that can be tuned separately, allowing the observation 
of two frequencies simultaneously. The 'backend' is a digi- 
tal autocorrelator, providing two simultaneous outputs with 
available band widths ranging from 4 to 256 MHz, which can 
be split into a maximum of 4096 channels. For the duration 
of this survey, the only observing method available was posi- 
tion switching: i.e. the observations were made by switching 
between the on-source position and a reference position with 
insignificant emission. 

Both the antenna and the digital correlator have fea- 
tures in common to those used in the Australia Telescope 
Compact Array, and further details can be found in the 
Special Edition of the Journal for Electrical and Elec- 
tronic Engineering, Austraha, Vol 12, No 2 (1992). Tech- 
nical details specific to the Mopra antenna can be found at 
www.narrabri.atnf.csiro.au/mopra/mopragu.pdf. 

An important aspect of the study reported here (involv- 
ing observations collected between 1994 and 2000) was the 
formulation of an intensity calibration scheme for the Mo- 
pra telescope. This calibration scheme is discussed in the 
appendix. 



2.1.2 The SEST Telescope 

The SEST is situated at the European Southern Observa- 
tory at La Silla in Chile, about 600 km north of Santiago. 
The SEST has a Cassegrain design and a parabolic collect- 
ing surface with a diameter of 15 metres, giving a half power 
beam width ranging from 57 arcsec at 86 GHz to 15 arcsec 
at 346 GHz. 

The SEST is fitted with three SIS receivers that collec- 
tively cover the frequencies 78 - 117 GHz, 128 - 170 GHz, 
215 - 270 GHz and 320 - 363 GHz. A polarization splitter 
allows the observation of two different frequencies simultane- 



^ Resulting from funding contributions from the University of 
New South Wales, the solid surface was extended over the full 
diameter in 1999, extending the millimetre-wave collecting area 
to 22 metres. 



ously, although the combinations of frequencies are limited 
to 3- mm and 2- mm, or 3-mm and 1.6-mm. 

Three Acousto-Optical Spectrometers (AOS) provide 
backends to the receivers. The high resolution AOS has a 
maximum bandwidth of 86 MHz, which can be divided into 
up to 2000 channels. 

For more instrumentation details see Booth et al. (1989) 
or The SEST Handbook (1997). 



2.2 Selection of the Survey Sample 

The aim of the selection process was to select around 25 
southern molecular clouds with strong molecular line emis- 
sion and evidence of complex chemistry for detailed study, 
and involved a two-step process. Firstly, previous molecular- 
line surveys of the southern Galactic plane (Whiteoak 
& Gardner 1974, Whiteoak & Gardner 1978, Gardner & 
Whiteoak 1978, Batchelor et al. 1981, Dickinson et al. 1982, 
Gardner & Whiteoak 1984, Peters et al. 1986) were used to 
choose 50 candidate molecular clouds south of declination 
-30 degrees. Secondly, these candidates were observed with 
the Mopra radio telescope between September and Novem- 
ber 1994, in the 3-mm transitions of CS, HCN and HCO+ 
at frequencies of 98.0, 88.6 and 89.2 GHz respectively. These 
transitions were selected firstly because these molecules are 
generally abundant in molecular clouds and secondly be- 
cause they were within the limited range of frequencies avail- 
able on the Mopra telescope at that time. These observa- 
tions were used to locate the line emission peaks, and select 
24 of the brighter clouds for detailed study. The sample in- 
cluded several (5) dark clouds and also some objects show- 
ing unusual intensity ratios of the three transitions. Note 
that these observations were made during commissioning of 
the receiver, and before the development of better positional 
calibration and intensity calibration. Consequently, later ob- 
servations of these transitions were obtained for use in the 
quantitative study (see table 3). 

The list of 24 Galactic molecular clouds selected is given 
in table 1, with the observed positions in both B1950 and 
J2000 coordinates. The two molecular clouds Orion MC and 
M 17, which are near the celestial equator and well stud- 
ied by radio telescopes in the northern hemisphere, have 
been included for comparison and calibration. Four sepa- 
rate pointing positions were chosen within the complex NGC 
6334, and two positions within G291. 3-0.7 (NGC 3576, RCW 
57). Therefore, including the two reference clouds (Orion 
and M 17), there are 30 target positions in table 1. 

The exact positions were chosen based on ^^CO 1-0 ob- 
servations at 110 GHz with the Mopra telescope in Septem- 
ber 1995. The clouds were mapped over a region of a few 
square arcminutes to determine the position of the peak, 
and the angular size of the ^^CO emission. The beam size 
was 37 arcsec, and the pointing errors 10-15 arcsec, with oc- 
casional discrepancies of up to 20 arcsec; the positions listed 
in table 1 are believed to be within 20 arcsec of the "CO 
condensation peak. Later observations centred on these po- 
sitions (table 3) have pointing errors better than 10 arcsec. 
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Table 1. Galactic molecular cloud positions observed. 



General Other Right Asc. Decl. Right Asc. Dccl. 

Designation Designation (J2000) {J2000) (B1950) (B1950) 



Orion M.C. 


Orion KL 


05 


35 


14, 


,5 


-05 


22 


29 


05 


32 


47. 


,0 


-05 


24 


22 


HH46 D.C. 




08 


25 


37, 


.8 


-51 


04 


00 


08 


24 


10, 


,6 


-50 


54 


08 


G265.1+1.5 


ROW 36 


08 


59 


26, 


,4 


-43 


45 


08 


08 


57 


38, 


,0 


-43 


33 


24 


G268.4-0.8 




09 


01 


51, 


,6 


-47 


44 


07 


09 


00 


09, 


,4 


-47 


32 


15 


Cham D.C. 




11 


06 


32, 


.6 


-77 


23 


35 


11 


05 


08, 


1 


-77 


07 


20 


G291. 3-0.7 


ROW 57 (E) 


11 


12 


04, 


.1 


-61 


18 


29 


11 


09 


56, 


1 


-61 


02 


09 


G291. 3-0.7 


ROW 57 (W) 


11 


11 


47, 


,4 


-61 


19 


28 


11 


09 


39, 


,6 


-61 


03 


09 


Coalsack D.C. 




12 


01 


27, 


,7 


-65 


07 


54 


11 


58 


53, 


,8 


-64 


51 


11 


G301.0+1.2 


RCW 65 


12 


34 


57, 


.8 


-61 


39 


56 


12 


32 


06, 


1 


-61 


23 


24 


G305.4+0.2 




13 


12 


33, 


.8 


-62 


33 


49 


13 


09 


21, 


,0 


-62 


17 


54 


G311. 6+0.3 




14 


04 


56, 


.0 


-60 


20 


21 


14 


01 


20, 


,0 


-61 


06 


00 


G322.2+0.6 


ROW 92 


15 


18 


38, 


,0 


-56 


39 


00 


15 


14 


47, 


,6 


-56 


28 


05 


Lupus D.C. 




15 


43 


02, 


,2 


-34 


09 


06 


15 


39 


51, 


,2 


-33 


59 


36 


G326.7+0.6 




15 


44 


44, 


.8 


-54 


06 


25 


15 


40 


55, 


,0 


-53 


57 


00 


G327.3-0.5 




15 


53 


06, 


.1 


-54 


35 


31 


15 


49 


13, 


,0 


-54 


26 


37 


G331. 5-0.1 




16 


12 


12, 


.3 


-51 


29 


02 


16 


08 


24, 


,2 


-51 


21 


20 


G333.0-0.6 




16 


21 


05, 


.5 


-50 


35 


25 


16 


17 


18, 


,3 


-50 


28 


18 


G333.4-0.4 




16 


21 


32, 


.9 


-50 


26 


32 


16 


17 


46, 


,0 


-50 


19 


27 


G333.6-0.2 




16 


22 


12, 


.1 


-50 


05 


56 


16 


18 


26, 


,0 


-49 


58 


54 


G345.5+1.5 




16 


59 


41, 


.9 


-40 


03 


13 


16 


56 


14, 


1 


-39 


58 


45 


G345.5+0.3 




17 


04 


28, 


,5 


-40 


46 


02 


17 


00 


59, 


,0 


-40 


41 


54 


NGC 6334 


NGC 6334 (S) 


17 


19 


55, 


.9 


-35 


57 


53 


17 


16 


34, 


,5 


-35 


54 


51 


NGC 6334 


NGC 6334 (CO) 


17 


20 


23, 


.4 


-35 


55 


00 


17 


17 


02, 


,0 


-35 


52 


00 


NGC 6334 


NGC 6334 (N) 


17 


20 


48, 


.5 


-35 


46 


33 


17 


17 


27, 


,4 


-35 


43 


35 


NGC 6334 


NGC 6334 (Nl) 


17 


20 


53, 


.4 


-35 


45 


32 


17 


17 


32, 


,3 


-35 


42 


35 


G348.7-1.0 


RCW 122 


17 


20 


05, 


.2 


-38 


57 


22 


17 


16 


38, 


,3 


-38 


54 


21 


G351.6-1.3 




17 


29 


12, 


.8 


-36 


40 


12 


17 


25 


49, 


,8 


-36 


37 


50 


G353.4-0.4 




17 


30 


24, 


.6 


-34 


41 


39 


17 


27 


05, 


,0 


-34 


39 


23 


Cor Aust D.C. 


R CrA 


19 


10 


18, 


,2 


-37 


08 


58 


19 


06 


56, 


,0 


-37 


13 


54 


M 17 


M 17 SW 


18 


20 


23, 


,1 


-16 


11 


36 


18 


17 


30, 


,0 


-16 


12 


59 



2.3 Survey Observations 

2.3.1 Mopra Observations 

Molecular transitions between 85 and 116 GHz were ob- 
served with the Mopra telescope between September 1994 
and September 1996 with the 15-metre solid millimetre sur- 
face, and in March and September 2000, after the resurfacing 
to 22 metres. Emission from 11 different molecules was ob- 
served (CO, CH3OH, C2H, CS, HCN, HC3N, HCO+, HNC, 
HNCO, OCS and SO), generally in more than one transition. 
Observations of more rare, isotopically substituted species 
of the molecules CO, CS, HCN, HCO+ and HNC were also 
obtained (^^C for ^^C and ^*0 for ^''O). The actual transi- 
tions observed are listed in table 2, together with the date 
of observation, the frequency and the telescope used. 

The single-sideband receiver temperature varied be- 
tween 85 and 130 K. All observations were obtained with 
a bandwidth of 64 MHz divided into 1024 channels, giving 
a velocity resolution between 0.44 and 0.32 km s~^ after 
Banning smoothing. 

Observations were made in position switching mode, 
generally with a separation of 30 arcmin; however, for some 
extended sources such as the dark clouds, 60-arcmin sep- 
arations were used. Periodic observations of SiO masers of 
known position were used to correct the pointing model, giv- 
ing a pointing accuracy of better than 10 arcsec for observa- 



tion obtained in 1995 or later. This includes all observations 
reported in table 3. 



2.3.2 SEST Observations 

The SEST was used in June 1999 to observe transitions be- 
tween 84 and 147 GHz. The system temperatures were 170 
- 250 K at 3 mm and 230-400 K at 2 mm wavelength bands. 
The beamwidth varied from 59 arcsec at 85 GHz to 35 arc- 
sec at 147 GHz. Note that although the effective diameter 
of SEST was the same as Mopra (15 m) the Mopra dish was 
'shaped' with different illumination resulting in a smaller 
beamwidth for a given frequency. A bandwidth of 42 MHz 
was used, divided into 1000 channels, giving velocity resolu- 
tion between 0.30 and 0.17 km s~^ after Hanning smoothing. 

The observations were made using two different ob- 
serving modes. Generally, dual-beam switching mode was 
utilised with a throw of 11 arcmin between the on-source and 
reference positions. However, for extended sources, where 
line emission was likely to be present in the reference posi- 
tion, frequency switching was used, with a frequency throw 
of 4 MHz. Periodic observations of SiO masers were used to 
check focus and pointing, giving a pointing accuracy better 
than 5 arcsec. Intensity calibration of the observations was 
performed on-line by the chopper wheel method (see Ap- 
pendix), with the spectra obtained being on the Ta* scale. 
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Table 2. Molecular transitions observed 



JVlolecule 


Xransition 




rvest Frequency 


Telescope 


uare 
Observed 


CO 


1-0 




115.27120 


Mopra 


1996 May 


"CO 


1-0 




110.20135 


Mopra 


1996 May 


ClSQ 


1-0 




109.78216 


Mopra 


1996 May 


CH30H 


5(-l)-4(0) E 




84.52121 


SEST 


1999 Jun 


CH30H 


15(3)-14(4) A- 




88.94009 


SEST 


1999 Jun 


CH30H 


2(-l)-l(-l) E 




96.73939 


Mopra 


1995 Nov 


CH30H 


2(0)-l(0) A+ 




96.74142 


Mopra 


1995 Nov 


CH30H 


2(0)-l(0) E 




96.74458 


Mopra 


1995 Nov 


CH30H 


2(1)-1(1) E 




96.75551 


Mopra 


1995 Nov 


CH30H 


3(l)-4(0) A+ 




107.01385 


SEST 


1999 Jun 


CH30H 


3(0)-2(0) E 




145.09375 


SEST 


1999 Jun 


CH30H 


3(-l)-2(-l) E 




145.09747 


SEST 


1999 Jun 


CH30H 


3(0)-2(0) A+ 




145.10323 


SEST 


1999 Jun 


C2H 


1-0 3/2-1/2 F= 


-l-l 


87.28415 


Mopra 


1995 Nov 


C2H 


1-0 3/2-1/2 F= 


:2-l 


87.31692 


Mopra 


1995 Nov 


C2H 


1-0 3/2-1/2 F= 


-1-0 


87.32862 


Mopra 


1995 Nov 


cs 


2-1 




97.98096 


Mopra 


1995 Nov 










SEST 


1999 Jun 


"cs 


2-1 




92.49429 


Mopra 22m 


2000 Mar 


cs 


3-2 




146.96904 


SEST 


1999 Jun 


HON 


1-0 F=l-1 




88.63041 


Mopra 


1995 Oct 


HON 


1-0 F=2-l 




88.63184 


Mopra 


1995 Oct 


HON 


1-0 F=0-1 




88.63393 


Mopra 


1995 Oct 


HC3N 


10-9 




90.97899 


Mopra 


1995 Nov 










SEST 


1999 Jun 


HC3N 


11-10 




100.07638 


SEST 


1999 Jun 


HC3N 


12-11 




109.17363 


Mopra 


1996 May 


HC3N 


15-14 




136.46440 


SEST 


1999 Jun 


HC3N 


16-15 




145.56094 


SEST 


1999 Jun 


HCO+ 


1-0 




89.18851 


Mopra 


1995 Oct 


H"CO+ 


1-0 




86.75429 


Mopra 22m 


2000 Oct 


HC"0+ 


1-0 




85.16225 


Mopra 


1996 Sep 


HNC 


1-0 F=0-1 




90.66345 


Mopra 


1995 Nov 


HNC 


1-0 F=2-l 




90.66357 


Mopra 


1995 Nov 


HNC 


1-0 F=l-1 




90.66365 


Mopra 


1995 Nov 


HN"C 


1-0 




87.09086 


Mopra 22m 


2000 Oct 


HNCO 


4(2,3)-3(2,2) 




87.89841 


Mopra 


1996 Sep 


HNCO 


4(2,2)-3{2,l) 




87.89862 


Mopra 


1996 Sep 


HNCO 


4(0,4)-3{0,3) 




87.92523 


Mopra 


1996 Sep 


OCS 


7-6 




85.13910 


Mopra 


1996 Sep 


OCS 


8-7 




97.30120 


Mopra 


1996 Sep 


OCS 


9-8 




109.46306 


Mopra 


1996 Sep 


OCS 


12-11 




145.94681 


SEST 


1999 Jun 


SO 


3(2)-2(l) 




99.29987 


Mopra 


1995 Nov 



2.4 Data Reduction 

All spectra were reduced using the CLASS (Continuum and 
Line Analysis Single-dish Software) package of the GILDAS 
(Crenoble Image and Line Data Analysis Software) working 
group. 

The spectra were Hanning-smoothed to increase the 
signal-to-noise ratio before baseline subtraction and Gaus- 
sian fitting to features. The spectral profiles generally had 
complex shapes, so one or more Gaussian components were 
fitted. 

For Mopra spectra, the peak and integrated intensities 
of the Gaussian components were then corrected for atmo- 
spheric attenuation and scaled to the SEST corrected an- 
tenna temperature (T^) scale (see Appendix for more details 



on intensity calibration of Mopra). The SEST corrected an- 
tenna temperatures (TJ) were then converted to main-beam 
brightness temperatures (Tmb) using the SEST main beam 
efficiency {B^f f / F^f f) which varied from 0.75 to 0.66 be- 
tween 87 and 147 GHz (Booth et al, 1989, SEST Observers 
Handbook). 



3 RESULTS 

Table 3 lists the observed molecular transitions, and the 
results of the Gaussian-component fits for each molecular 
cloud. An excerpt from table 3, showing results for the 
molecular cloud G265.5-I-1.5, is given in the body of this 
paper. The complete table can be found in the electronic 
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version of the paper and on the VizieR data base of astro- 
nomical catalogues at the Centre de Donnees astronomiques 
de Strasbourg (CDS) website. The first two columns list 
the molecule and the transition observed. The next three 
columns are the Gaussian fit parameters (central velocity, 
full-width at half maximum and peak expressed as main 
beam brightness temperature Tmb)- Multiple Gaussians 
were used to fit complex profiles. However, interpreting com- 
ponents in terms of separate clouds is difficult because in 
some cases, double-peaked profiles fitted as two Gaussians 
may be a broad emission line with a self-absorption dip. The 
velocities of the HCN hyperfine components are for the J = 
2-1 component with the F — 1-0 and F = 1-1 line centres 
fitted using the same velocity scale. These latter hyperfine 
transitions are shown in brackets in the table. (This is be- 
cause the components are often blended, so it was not feasi- 
ble to velocity correct the hyperfine components separately.) 
The 6th column shows the integrated line emission (/ Tmb 
dv) for individual Gaussian components, while the 7th col- 
umn shows the total line emission in the transition, obtained 
by adding the individual Gaussian components; in the case 
of HCN this total includes the emission in the three hyper- 
fine components. The main beam brightness temperature 
(Tms), is the true brightness temperature of the source, if 
the source just fills the main beam. If the source angular size 
is known, then this temperature can be deconvolved with the 
source size to give the actual source brightness temperature. 

The 8th column gives the uncertainty in the integrated 
brightness temperature (column 7). For Mopra observations, 
the scatter of repeated observations of a source was used to 
make an estimate of the accuracy of the integrated bright- 
ness temperature of the lines. The error in integrated bright- 
ness was considered to be made up of two parts adding in 
quadrature - a constant percentage error B, which domi- 
nates at high brightness, and a constant absolute error A, 
due to noise in the spectra, which dominates at low bright- 
ness. From the scatter of repeated observations of strong 
lines {JTmb dv > 5 K km s~^) we determine B = 12 %, 
and from repeated observations of weak lines {J Tmb dv < 
5 K km s~^) we determine A = 0.24 K km s~^. For a typical 
linewidth of 4 km s~^, this value of A corresponds to peak 
brightness 0.06 K and is consistent with the contribution to 
line fitting from the noise in the spectra. 

The calibration uncertainties of the SEST observations 
are assumed to be 10 % at 3 mm and 15 % at 2 mm (SEST 
Observers Manual 2003; Heikkila 1998), and the total uncer- 
tainty the sum in quadrature of the calibration uncertainty 
and the formal errors in fitting the lines using the CLASS fit- 
ting routine (GAUSS). The Mopra and SEST uncertainties 
turn out to be similar (listed in column 8 of table 3). 



4 DISCUSSION 

Some general properties of the line emission and profiles are 
discussed in this section, including the ratios of the emission 
from different isotopomers of several molecules. In addition, 
the line intensities have been used to calculate molecular 
abundances using both LTE and LVG models, which we 
present in a separate paper (Hunt Cunningham et al. 2004) . 



4.1 Spectral Line Intensities and Profiles 

Figure 1 shows CS 2-1 spectra towards each position listed 
in table 1. Each spectrum is plotted over a velocity inter- 
val of 30 km s"'^. This transition was chosen to illustrate 
the characteristic line profile in each molecular cloud as it 
is ubiquitous, generally has good signal to noise, and tends 
to have a 'typical line profile. On the other hand, the CO 
spectra for each molecular cloud often show emission fea- 
tures covering a significantly wider range of velocities than 
found for other molecules, as can be seen in figure 2, which 
shows spectra for a number of different transitions towards 
the molecular cloud G265.1-I-1.5 (RCW 36). The spectra to- 
wards this molecular cloud illustrate features common to 
many of the molecular clouds. 

Self-absorption is evident in the ^'^CO spectra towards 
about a third of the molecular clouds, and can also be seen 
in some spectra of '^^CO and C^*0. Transitions that are 
believed to be self-absorbed are marked in table 3 by an 
asterisk. Self-absorption is also evident in some clouds in 
transitions of CS, HCN, HCO"'' and SO, but not gener- 
ally in transitions of OCS, HC3N and CH3OH, implying 
that these latter transitions are more likely to be optically 
thin. In addition to the presence of self-absorption, the spec- 
tra, particularly those of ^^CO, often show a complex emis- 
sion/velocity structure that cannot be adequately fitted by a 
single Gaussian. As expected, emission from the dark clouds 
in the sample (HH 46, the Coal Sack, the Chameleon, Lu- 
pus and Corona Australis) is distinguished by narrow line 
widths compared to other clouds in the sample (see figure 
!)• 

High-velocity wings and asymmetric line profiles are 
seen in the spectra of a number of molecular clouds (e.g 
spectra of CO, CS, HCN and SO observed towards the po- 
sitions NGC 6334(CO) and NGC 6334 (N)) suggesting that 
outfiows, presumably from young stellar objects (YSOs) , are 
present. Figure 3 shows some examples of spectra with line 
wings and asymmetries. 



4.2 Isotopomer Intensity Ratios 

The ^■^CO emission is found to be optically thick for all 
molecular clouds in this sample. Table 4 lists the emission 
intensity ratios for ^^C'}^C isotopomers of CO, CS, HCN, 
HNC and HCO^, as well as the sample standard devia- 
tion (cr„_i), the minimum and maximum values of the ra- 
tio, and the number of molecular clouds for which the ra- 
tio has been calculated. The terrestrial abundance ratio of 
the atoms ^^0:^*^0 is 89, the ratio of the abundances of the 
molecules ^'^COi^^CO in nearby molecular clouds is found 
to be around 50; however, significant variations in this ra- 
tio have been found in individual molecular clouds, and it 
is likely that the ratio varies with galactocentric distance. 
Langer and Penzias (1990) have found a value of 57 for 
molecular clouds at the Sun's distance from the Galactic 
Centre, with values ranging from 25 near the centre of the 
Galaxy to 70 at a distance of 12 kpc. In contrast, the mean 
ratio for the intensity of ^^CO emission (hereafter I^'^) to 
emission from ^^CO (hereafter l") in table 4 is 3.1 ± 0.2 
for this sample of southern molecular clouds, implying that 
the i^CO emission is optically thick for all molecular clouds 
of the sample. This is very close to the value for I^^: I^^ 
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Table 3. Fitted line intensities from Mopra and SEST observations for the molecular cloud G265.5+1.5 (RCW 36). The remainder of 
the table appears in the electronic version of this paper. 



Gaussian Components of fit Total Emission 

Molecule Transition Centre dv Ta/s A/sdv (5 J Tji/Bdv 

(km s-i) (km s-i) (K) (K km s"!) (K km s"!) (K km s'^) 



12CO 


1-0 


1.4 


1.3 


14.25 


19.73 


19.7 


2.3 


"CO 


1-0 


4.0 


0.9 


6.90 


6.62 


6.6 


0.8 


c^°o 


1-0 


4.5 


0.6 


1.08 


0.69 


0.7 


0.3 


C2H 


1-0 3/2-1/2 F=l-0 


4.4 


2.0 


0.97 


0.25 


0.7 


0.3 




1-0 3/2-1/2 F=2-l 


1.1 


0.8 


0.39 


0.43 






OS 


2-1 


5.3 


0.9 


1.04 


0.99 


1.0 


0.3 


CS 


3-2 


5.2 


0.7 


0.56 


0.45 


0.4 


0.2 


HCN 


1-0, F=2-l 


5.3 


1.1 


0.87 


1.02 


1.8 


0.3 




1-0, F=l-1 


(10.1) 


1.0 


0.38 


0.40 








1-0, F=0-1 


(-1.7) 


1.0 


0.33 


0.35 






Ri^CN 


1-0 




2.4 


0.31 


0.78 


0.8 


0.2 


HNC 


1-0 


5.2 


0.9 


0.72 


0.67 


0.7 


0.3 


HN"C 


1-0 






- 




<0.2 




HNCO 


4(0,4)-3(0,3) 










<0.10 




HCO+ 


1-0 


5.3 


1.1 


1.40 


1.63 


1.6 


0.3 


Hi3CO+ 


1-0 


4.9 


1.0 


0.39 


0.41 


0.4 


0.1 


HCi*0+ 


1-0 










<0.14 




HC3N 


10-9 






_ 




<0.04 




HC3N 


11-10 






- 




<0.07 




HC3N 


12-11 










<0.06 




HC3N 


15-14 






- 




<0.08 




HC3N 


16-15 










<0.08 




CCS 


7-6 










<0.14 




DCS 


8-7 










<0.09 




OCS 


9-8 










<0.23 




DCS 


12-11 










<0.09 




CH3OH 


2(0)-l(0) A-f 










<0.31 




CH3OH 


3(-l)-2(-l)E 










<0.13 




CH3OH 


15(3)-14(4)A- 










<0.04 




SO 


3(2)-2(l) 


5.2 


0.7 


1.29 


0.99 


1.0 


0.3 



of 3.0 ± 0.9 found for a sample of northern Galactic-plane 
GMCs between 1 = 34° and 1 = 51° (Polk et al. 1988) sug- 
gesting that there is no systematic difference in this ratio 
between the northern and southern Galactic plane. Doty & 
Neufeld (1997) have modelled I^^: I^^ for the dense quiescent 
star-forming cores of molecular clouds, assuming an approx- 
imately terrestrial abundance ratio of 100:1 for ^^CO:^''CO. 
They have predicted an I^^: I^^ ratio of between 2 and 6, 
which is consistent with the findings of this study, and again 
suggests that the ^^CO emission has high optical depth for 
all clouds in this sample. 

Table 5 lists the same values for the intensity ratio of 
^•^C: isotopomers. The mean ratio of I^'^ to emission 
from C^*0 (hereafter l") is 6.9 ± 0.4. This is somewhat 
higher than the value for this ratio of 5.5 that would be 
expected for optically thin emission if terrestrial ratios for 
atomic ^2C:"C and ^'^O: ^**0 were assumed. The values for 



1:1 for individual clouds range from 2.5 to 15.2, so in 
some molecular clouds the ratio is far higher than expected 
from the terrestrial ratio of the isotopes. 

The ratio I^^: I^* found in this study is 25 ± 3, with a 
range of 4.2 to 98. Doty & Neufield (1997), assuming ter- 
restrial ratios, predicted that I^^: I^* should be between 5 
and 20, assuming that ^^CO and C^*0 were coextensive. A 
possible reason for the high ratios found in this study is that 
the emission from ^^CO arises from a larger volume of gas 
than does emission from C^*0. This suggestion is supported 
by the average line-widths for ^'^CO and C^*0 emission for 
the sample. The average line-width for ^^CO is 8.7 km s~^, 
while the average line-width for C^^^O is 4.1 km s~^, sug- 
gesting that ^■^CO emission arises from a larger volume of 
the gas (assuming that line-width is proportional to radius, 
as if will be if there is a linear velocity gradient across the 
molecular cloud). However, note that this difference in av- 
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Figure 1. CS 2-1 spectra towards each of the molecular clouds in table 1. This transition shows features characteristic of the many of 
the other molecular transitions observed. The name of each molecular cloud is shown above each spectrum. 



Table 4. Emission intensity ratios for ^'^C:^^C from isotopomers of CO, CS, HCN, HNC and HCO+. 





CO 


CS 


HCN 


HNC 


HCO+ 






H^'-'CN 


HN^'-'C 


m-^co+ 


Mean 


3.1 


12.0 


7.8 


11.9 


6.8 


Standard dev. 


1.9 


5.7 


5.4 


5.6 


7.6 


Minimum 


1.7 


6.1 


0.9 


2.7 


0.3 


Maximum 


10.0 


30.2 


25.1 


21.7 


33.3 


Number of Obs. 


33 


19 


25 


20 


28 



erage line widths may be due to optical depth effects. If the 
^^CO emission is optically thin in the wings of the line pro- 
files then the C^*0 emission in the line wings will not be 
detected, assmning a terrestrial abundance of ^^CO:C^*0 of 
490:1. 

The average line-width of ^"^CO is 5.3 km s~^, simi- 
lar to that for C^^O, suggesting that the C^**0 and ^^CO 



emission arise from a similar volume of gas. If this is the 
case, then the optical depth of the ^'^CO emission is likely 
to be even higher than that implied by the emission ratios 
in Table 4. The critical densities for ^^CO and C^**0 are 
lower than that for ^^CO, implying that the difference in 
volume of the emitting gas is due to a lower abundance of 
the more rare isotopomers in outer regions of the molecu- 
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Figure 2. Spectra towards the molecular cloud G265.5+1.5. The transition is shown in the top left-hand corner of each box. The 
differences between the line profiles of the different transitions are typical of many of the molecular clouds in the sample. 



lar cloud rather than to varying excitation conditions. The 
unexpectedly high values of I^'^iP* can be understood if the 
^^CO and C^*0 emission are not completely co-extensive, 
with the C^*0 emission arising from a smaller volume of 
gas. The relative average line widths (4.1 km s~^ for C^*0, 
5.3 km s~^ for ^^CO) provide some support for this asser- 
tion. It is possible that self-shielding (from dissociating UV 
radiation both within and incident on the molecular cloud) 
in the more abundant isotopomers accounts for their appar- 
ently wider distribution within the molecular gas (Langer & 
Penzias 1990). 

The average ratio of emission from H^^CO'*' to that 
from HC^*0+ is 15.3 ± 3.5 (Table 5). Once again, this is far 
higher than would be expected from optically thin emission 



from both isotopomers if terrestrial ratios of atomic C: 
"C and ^^O: ^'^O apply. The high ratio suggests that the 
emission may arise from a smaller volume of gas 
than the R^^CO+ emission, once again possibly due to self- 
shielding in the more abundant isotopomer (H^"^CO^). The 
higher ratio for emission from H^^CO+:HC^*0+ compared 
with that from ^'^CO:C^*0 may also support this sugges- 
tion. HC^*0+ can be expected to be far less abundant than 
C^*0, and consequently full self-shielding from dissociat- 
ing UV radiation will occur at a significantly smaller cloud 
radius. On the other hand, the high ratio may be due to 
chemical fractionation effects. Chemical fractionation (the 
preference for a molecule consisting of one isotope over an- 
other in a chemical reaction) can lead to different isotope 
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Figure 3. Examples of spectra showing line wings and asymmetric line profiles, suggesting the presence of outflow. 



Table 5. Emission intensity ratios for ^''CO:C^^O from iso- 
topomers of CO and HCO"*". 









Mean 


6.9 


15.3 


Standard dev. 


2.6 


13.7 


Minimum 


2.5 


6.3 


Maximum 


15.2 


60.0 


Number of Obs. 


33 


15 



ratios in the molecular products of a reaction with respect 
to the reactants. More complex products have a higher ra- 
tio of the lighter to the heavier isotope (e.g a higher ^^O: 
^**0 ratio) (Engel & Macko 1997), so HC^®0+ may be less 
abundant with respect to its ^'^C isotopomer than C^*0. 

The average emission ratios of the ^'^C to ^'^C iso- 
topomers of CS, HCN, HNC and HCO+ (Table 4) suggest 
that the ^^C isotopomer of these molecules is generally opti- 
cally thick in dense molecular clouds, although less so than 
^■^CO. A surprising result is that the minimum value for the 
emission ratio HCN:H^''CN is 0.9, and the minimum value 
for HCO^:HC130^ is 0.3. Both these low values occur in 
the same molecular cloud, the Chameleon dark cloud. No 
other molecular clouds have values less than unity for these 
ratios. In fact, all ratios for emission from ^^C isotopomers 
to ^^C isotopomers in the Chameleon dark cloud are low 



compared to other molecular clouds. The reason for this is 
not clear. It is possible that the Chameleon dark cloud has 
a low elemental ratio of ^^C to ^''C, but further observations 
of the Chameleon dark cloud will be required to clarify the 
situation. 



5 SUMMARY 

Transitions of 11 different molecules and several of their iso- 
topomers have been observed in a sample of 26 southern 
molecular clouds. The observational results are as follows: 

1. The CO spectra often show a significantly wider range 
of velocities than found for other molecules. 

2. Self-absorption is present in about a third of the 
CO profiles, and in some profiles of "CO, C^*0, CS, HCN, 
HCO+ and SO. 

3. Self-absorption is not generally observed in transi- 
tions of OCS, HC3N and CH3OH, suggesting that these 
transitions are optically thin in most molecular clouds of 
the sample, or that these transitions only exist in the higher- 
excitation, inner parts of the clouds. 

4. CO is optically thick in all molecular clouds of 
the sample. The average emission ratio for ^■^CO:^'^CO is 
3.1 ± 0.2, very close to values found in surveys of northern 
hemisphere molecular clouds, and predicted by models. 

5. The average emission ratio for ^^CO:C^*0 is 6.9±0.4, 
higher than would be expected for optically thin emission if 
Solar System ratios for ^^C:^*0 were assumed. This may be 
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because O occupies a smaller volume of gas than does 
"CO, due to self-shielding of "CO. 

6. The average emission ratio of H"C0+:HC"0+ is 
15.3 ± 3.5, also higher than would be expected for optically 
thin emission if Solar System ratios for "C:"0 were as- 
sumed. Self-shielding in H^'^CO^ is a possible explanation, 
however chemical fractionation effects may also be respon- 
sible for the high ratio. 

7. CS, HCN, HNC and HCO+ are generally found to be 
optically thick for the molecular clouds of the sample. 

8. The Chameleon molecular cloud has ^^C:^^C iso- 
topomer ratios of less than unity for IfCN and HCO^, and 
has low ^^C:^''C isotopomer ratios for all molecules. Further 
observations will be required to clarify the reasons for this 
situation. 



ACKNOWLEDGMENTS 

This research forms part of the PhD thesis of MRHC under- 
taken at the University of Western Sydney. MRHC acknowl- 
edges support from an Australian Postgraduate Award and 
thanks DIST Access to Major Research Facilities Scheme for 
travel support to SEST. MRHC also gratefully acknowledges 
support from the ATNF, including free accommodation and 
use of computing facilities for the duration of this work. The 
Mopra telescope is part of the Australia Telescope which is 
funded by the Commonwealth of Australia for operation as 
a National Facility managed by CSIRO. 



REFERENCES 

Batchelor, R.A., McCuUoch, M.G. & Whiteoak, J.B. 1981, MN- 
RAS, 194, 911 

Booth, R.S., Delgado, G., Hagstrom, M., Johansson, L.E.B., Mur- 
phy, D.C., Olberg, M., Whyborn, N.D., Greve, A., Hannson, 
B., Lindstrm, CO. & Rydberg, A. 1989, A&A 216, 315 

Dickinson, D.F., Gulkis, S., Klein, M.J., Kuiper, T.B.H., Batty, 
M., Gardner, P.P., Jauncoy, D.L. & Whiteoak, J.B. 1982, AJ, 
87, 1202 

Doty, S.D. & Neufeld, D.A. 1997, ApJ, 489, 122 
Engel, M.H., & Macko, S.A. 1997, Nature, 389, 265 
Gardner, P.P. & Whiteoak, J.B. 1978, MNRAS, 183, 711 
Gardner, P.P. & Whiteoak, J.B. 1984, MNRAS, 210, 23 
Heikkila, A. 1998, Technical Report No 343 (Thesis) Chalmers 

University of Technology 
Hunt Cunningham, M.R., Whiteoak, J.B., Jones, P.A., & White, 

G.L., 2004, MNRAS, in preparation 
Kutner, M.L. & Ulich, B.L. 1981, Ap. J., 250, 341 
Langer, W.D. & Penzias, A.A., 1990, Ap.J., 357, 477 
Peters, W.L. Porster, J.R. Gardner, P.P. Whiteoak, J.B. & 

Kuiper, T.B.H. 1986, PASA, 6, 504 
Polk, K.S., Knapp, G.R., Stark, A.A. & Wilson, R.W. 1988, Ap.J., 

332, 432. 

Rohlfs, K., & Wilson, T.L., 1996. Tools of Radio Astronomy, 2nd 

Edition, (Berlin: Springer- Verlag) 
SE ST Handbook, The. 1997, 

"http;// www. ls.eso.org/lasilla/Telescopes/SEST/SEST. html' 
van Dishoeck, E.P. & Blake, G.A. 1998, ARA&A, 36, 317 
Whiteoak, J.B. & Gardner, P.P. 1974, A&A, 37, 389 
Whiteoak, J.B. & Gardner, P.P. 1978, MNRAS, 185, 33 



APPENDIX: CALIBRATION OF THE MOPRA 
TELESCOPE 

The accurate intensity calibration of the observations re- 
ported here is important, as they have been used to cal- 
culate abundances (Hunt Cunningham, in preparation). At 
millimetre wavelengths there is significant emission and ab- 
sorption of radiation by the atmosphere, and the compo- 
nent due to water vapour in the atmosphere can vary over 
timescales as short as a few minutes. One way to correct 
for these optical depth fluctuations is the 'chopper wheel' 
method, where an ambient temperature blackbody source 
is periodically inserted into the antenna beam as a refer- 
ence (see e.g. Rohlfs & Wilson 1996; Kutner & Ulich 1981). 
This method as applied to both the Mopra telescope and 
the SEST is outlined here. More information on the terms 
and definitions used below can be found in Rohlfs & Wilson 
(1996). 

The antenna temperature directly measured at the tele- 
scope (Ta) needs to be corrected by various factors to ac- 
count for atmospheric attenuation, aperture efficiency, far- 
sidelobe contributions (back scattered radiation) and elec- 
tronic losses and noise (receiver temperature). A quantity 
known as the calibrated antenna temperature (Tji) is cor- 
rected for atmospheric attenuation, backscattered radiation 
and receiver temperature, while what is known as the cor- 
rected antenna temperature (TJ) is also corrected for the 
forward efficiency of the telescope (Feff), so — T^/Feff. 

The forward beam includes the sidelobes in the forward 
direction, so a better parameter to use is the main beam 
brightness temperature (Tmb), where 

Tmb = T^/Beff = T'XFeff/Beff. 

The SEST observations were made with on-line chop- 
per wheel calibration system, and the observations were cor- 
rected by the on-line software to the corrected antenna tem- 
perature TX, so only the (frequency dependent) correction 
Feff/Beff needed to be applied to get Tmb- 

The Mopra system did not have such on-line chopper 
wheel or software, so a 'blackbody paddle' was inserted every 
few minutes to measure the gated total power (GTP) which 
could then be used off-line to calibrate the antenna tem- 
perature Ta for atmospheric attenuation and instrumental 
parameters. 

Let Th be the receiver temperature, Ts the source tem- 
perature, Ta the ambient temperature of the atmosphere 
and T the optical depth of the line of sight through the at- 
mosphere (r = To sec{Z), where Z is the zenith angle and ro 
the optical depth vertically through the atmosphere). The 
on-source signal 5* is then 

S = C[Tr + Tail - exp(-r)) + Ts exp(-r)], 

where C is a constant, and the off-source reference R is 

R = C[TR + Ta{l-exp{^r))]. 

The usual position switching, subtracting the reference 
from the signal, gives 
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S CTs exp(-r) 

correcting for the receiver temperature and the atmosphere 
emission, but not for the atmospheric absorption of the 
source emission. 

The paddle reference can be considered as a black body, 
at the ambient temperature Ta- The power Pr at the off- 
source position is 

Pr = C'[Tr + Tail - exp(-T))] 

where C' is another constant, and the power Pp with the 
black-body paddle in the telescope beam is 

Pp = C'[Tii + Ta] so Pp-Pr = C'Ta exp(-r). 

This can be used to correct the observations for atmo- 
spheric attenuation by 

{S - R)/{Pp - Pr) = {C/C')iTs/Ta) or 

Ts = {C'/C)Ta{S - R)/{Pp - Pr) 

where C and C' are constants of the telescope hardware and 
software to be determined and calibrated out. 

A complication of the Mopra spectral calibration here 
is that the observing program SPECTRA did not output 
the difference S — R oi the signal and reference spectra, but 
the quotient 

Q = K[{S/R) - 1]Tr = KTr{S - R)/R 

where K is a, constant required from the instrumental cali- 
bration. 

The standard calibration method equation is modified 

to 

QPr/{Pp - Pr) = KTr{S - R)Pr/[R{Pp " Pr)] 

= KTRTs/Ta or 

Ts = Q{Ta/KTR)Pr/{Pp ~ Pr). 

The Mopra telescope was not calibrated absolutely, for 
the period that these observations were made, so the factor 
K was significant, and we had to do a relative calibration by 
comparing Mopra and SEST observations of the calibrator 
sources Orion MC and M17. The integrated line intensities 
(in K km s~^) were used for 43 lines between 86 and 115 GHz 
observed at Mopra several times over a range of weather 
conditions. Before calibration, the scatter of repeated obser- 
vations in Q about the mean for that line was rms 17 % 
with some data points more than 30 % from the mean, in- 
dicating the effect of variable atmospheric absorption. After 
calibration, the scatter about the mean was 7 %, showing 
the effectiveness of the correction for absorption. Figure 4 
shows the uncalibrated Mopra integrated antenna tempera- 
tures [J Qdv) plotted with the SEST integrated, corrected 
antenna temperatures (jTXdv), while figure 5 shows the 
same data after the calibration procedure outlined above. 

The Mopra spectra quotient Q values were converted 
to the SEST corrected antenna temperature scale with 
a factor 1.62 ± 0.06 accounting for the instrumental factors 
in the Mopra SPECTRA Q scale. 

This paper has been typeset from a Tj^X/ ffl^jX file prepared 
by the author. 



Uncorrected Mopra vs SEST Intensities 
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Figure 4. Uncorrected Mopra integrated brightness tempera- 
tures (Q values) plotted with SEST calibrated integrated bright- 
ness temperatures. The uncalibrated Mopra values show a large 
scatter of around 40 %, due to uncalibrated atmospheric attenu- 
ation, as well as a scale error due to unknown instrumental pa- 
rameters when compared to the calibrated SEST data. 



Corrected Mopra vs SEST Intensities 
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Figure 5. The corrected Mopra TJ values, show a much tighter 
correlation with the SEST values (c.f. figure 4), with scatter 
around 20 %, and the comparison enables the scale factor of 1.62 
to be determined. 
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